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We study slow-subdiffusion in comparison to subdiffusion. Both of the 
processes are treated as random walks and can be described within con- 
tinuous time random walk formalism. However, the probability density of 
the waiting time of a random walker to take its next step ui(t) is assumed 
over a long time limit in the form ui(t) ~ l/t a+1 for subdiffusion, and 
in the form w(t) ~ h(t)/t for slow-subdiffusion [h(t) is a slowly varying 
function] . We show that Green functions for slow-subdiffusion and subdif- 
fusion can be very similar when the subdiffusion coefficient D a depends on 
time. This creates the possibility of describing slow-subdiffusion by means 
of subdiffusion with a small value of the subdiffusion parameter a, and the 
subdiffusion coefficient D a varying over time. 

PACS numbers: 05.40.-a, 05.60.-k, 05.40.Fb 



1. Introduction 

Anomalous diffusion is usually described within a framework of the con- 
tinuous time random walk in which the probability densiy of finding a ran- 
dom walker at time t at point x, P(x,t), depends on the probability den- 
sity of the waiting time of the random walker to take its next step cj(t) 
and the probability density of jump length \(x). When u:(t) and A(x) 
are independent then the probability density P(x, t) reads in terms of the 
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Laplace transform C{uj(t)} = d)(s) = / °° e st f(t)dt and the Fourier trans- 
form T{X(x)} = X(k) = e ikx f{x)dx [I] 

P(k,s)= 1 -^± L_. (1 ) 

s l-u(s)\(k) 

As far as we know, an inverse transform of the above equation in the most 
general case has not been found yet, with the expection of a few very special 
cases. For this reason, Eq. ([1]) is usually considered within the limit of small 
values of s and k. For subdiffusion, it is assumed that the first moment of 
(j(t) (the average value) is infinite and that the Laplace transfom w(s) has 
the following form for small values of s 

u{s) 1 - r a s a , (2) 

where subdiffusion parameter a obeys < a < 1 and r a is a positive 
parameter, whereas all the moments of the natural order of A(x) are finite. 
Assuming that X(t) is a symmetric function, the Fourier transform X(k) has 
the following form for small values of k 

X(k) = 1 - ^ 2 y , (3) 

where a 2 is the second moment of X(x). Computing the inverse Laplace and 
Fourier transforms leads to co(t) being proportional to l/t 1+a (therefore 
uj(t) is refered to as a heavy-tailed distribution) and X(x) in the form of 
Gauss distribution. A subdiffusion equation which is obtained as a result 
of continuous time random walk formalism is a linear differential equation 
with a fractional time derivative. 

A random walk in which u(t) is a superheavy-tailed distribution has 
been recently studied [21 [3] . The superheavy-tailed distribution means that 

u(t) ~ M , ( 4 ) 

where h(t) is a slowly varying function, i.e. h(xt)/h(t) —> 1 when t —> +oo 
for all positive values of x- All the fractional moments of ui(t) ({oj c (t)) = 
J x c uj(t)dt, c > 0) are infinite. Such a process qualitatively differs from 
subdiffusion and is usually called slow-sub diffusion. It is rather difficult to 
find an equation with fractional derivatives which describes such a process 
since an order of the derivative would be very close to zero. 

In this paper we will compare Green functions for subdiffusion and slow- 
subdiffusion and we will also show that slow-subdiffusion could be described 
by subdiffusion with small values of subdiffusion parameter a and with 
subdiffusion coefficient D a which depends on time. 
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2. Green functions 

Green functions for subdiffusive transport in an unrestricted system is 
well known. They can be obtained after substituting Eqs. ([2]) and ([3]) into 
Eq. ([1]) and then computing the inverse Laplace and Fourier transforms. 
After calculation we obtain 

P ( X ^) = TT-L=/a/2-l, a /2 (^-7^=) > ( 5 ) 



where [1] 

1 ~ 1 / a 



M*> «) = 7T^ E fcIr( _ fc/ 3 _ „) ("^ J • ( 6 ) 

The Green function for slow-sub diffusive transport in an unrestriced 
system was derived in [2] and reads over a long time limit 



P s (x,t) = ^e-WVm^ , (7 ) 

where V(t) is the complementary cumulative distribution function of waiting 
times V(t) = 1 - /q w(i)dt = J t °° uj(t)dt, V(0) = 1 and V(+oo) = 0. In the 
further consideration we take cj(t) in the following form 

_ (r - 1) In"' 1 r, 
{ ) ~ {r, + t)]n*(T, + t) ' (8) 



where r > 1 and 77 > 1, which provides [2] 



V(t) 



In r/ 
ln(r/ + 1) 



-ir-1 



(9) 



3. Comparision between the Green functions for 
slow subdiffusion and subdiffusion 

In Figs. [IHS] we present the comparision between Green functions for 
slow-sub diffusion and subdiffusion calculates from the formulae given by 
Eqs. (J7|) and ([5|), respectively. We are interested in testing if we can choose 
the parameters in such a way that functions (7) and (8) cover function (5). 
We arbitrarily chose the values of the parameters r, a 2 , r\ and a. In order to 
match Green functions for slow-sub diffusion and subdiffusion for all times, 
the subdiffusion coefficient D a has to depend on time. Therefore, subdiffu- 
sion coefficient D a was treated as a fitting parameter. All calculations were 
performed for r = 2.34, a 2 = 0.1 and a = 0.01 and the values of the rest 
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of the parameters were given in the legend and the figures' captions. It is 
interesting to find the dependance of the subdiffusion coefficient on the rest 
of parameters for the given value of t. The preliminary studies showed that 
the subdiffusion coefficent most strongly depends on rj. In Fig. [5]we present 
the dependence between the fitting parameter D a and the parameter rj. 




X 



Fig. 1. The comparision between Green functions for subdiffusion [Eq. (0] and 
slow-subdiffusion [Eq. The latter Green function was calculated for rj = 1.01. 
The values of time and subdiffusion coefficient are given in the legend. The lines 
represents the Green functions for slow-subdiffusion, the symbols — for subdiffu- 
sion. All quantities are given in arbitrarily chosen units. 



As can be seen in Figs. HHSl Green functions for subdiffusion and slow- 
subdiffusion are in quite good agreement for all values of rj for the subd- 
iffusion coefficient D a depending on time. Moreover, we can observe that 
the values of the subdiffusion coefficient depend strongly on the values of 
parameter rj. In Fig. [5] we can observe that for higher values of rj we have 
obtained lower values of D a . Therefore, we conclude that it is possible to 
describe slow-subdiffusion as subdiffusion with small values of subdiffusion 
parameter and a subdiffusion coefficient which depends on time. A more 
detailed description of the subject presented in this paper will be present 
elsewhere. 



lewandowska printed on January 22, 2013 



■5 




X 



Fig. 2. The same as in Fig. [TJ but for -q — 3.54. 

4. Final reamrks 

In this paper we have compared Green functions obtained for subdif- 
fusion and slow-subdiffusion. It seems that it is possible to descibe slow- 
subdiffusion by means of subdiffusion with a small value of subdiffusion pa- 
rameter and the subdiffusion coefficient depending on time although these 
two processes seem to be qualitatively different from each other. We have 
encountered a similar situation when we have approximated the solution to 
nonlinear sundiffusion equations by the solution to the subdiffusion equation 
with a fractional time derivative with the subdiffusion coefficient depending 
on time. 

It would be interesting to study the dependence of the subdiffusion co- 
efficient on time or the dependence of the agreement between the Green 
functions for slow-subdiffusion and subdiffusion on the subdiffusion param- 
eter a. Moreover, investigation like this could allow one to find a more 
transparent physical interpretation of slow-subdiffusion. 
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Fig. 4. The dependence of the subdiffusion coefficient D a on the parameter rj for 



